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Abstract 

We compute the eigenf unctions, energies and Betlie equations for a class of gener- 
alized integrable Hubbard models based on gl{n\m) © gl{2) superalgebras. The Bethe 
equations appear to be similar to the Hubbard model ones, up to a phase due to the in- 
tegration of a subset of 'simple' Bethe equations. We discuss relations with AdS/CFT 
correspondence, and with condensed matter physics. 
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1 Introduction 



The Hubbard model was first introduced as a model to describe the effects of the strongly 
correlated d-electrons in transition metals [1,2]. It was shown that this model possesses a 
Mott metal-insulator phase transition [3, 4] and is relevant to the studies of high-T^ super- 
conductivity [5,6]. Few exact results are known for the two or three-dimensional Hubbard 
model but nevertheless these models are still actively investigated. 

In contrast, the one-dimensional Hubbard model is integrable and was first diagonalized 
by means of coordinate Bethe Ansatz by Lieb and Wu [7] . However its integrable structure is 
rather complicated in comparison with the usual formalism of the spin chains (for a review 
see e.g. [8] and references therein). The Hubbard model i?-matrix was first introduced 
by Shastry [9] and Olmedilla et al [10] (by coupling two XX model /^-matrices with U- 
interaction term). The proof of the Yang-Baxter relation was given by Shiroishi and Wadati 
[11]. A lot of extended Hamiltonians were also proposed to connect the model with high-T^ 
superconductivity effects [12-14] (and references therein). 

Generalization of one-dimensional Hubbard model to gl{n) case via i?- matrix was made 
by Maassarani et al [15]. Extension of Maassarani's approach to superalgebras gl{n\m) 
was given in [16], and a "universal" presentation (including a generalization of Maassarani's 
approach) was proposed in [17]. It is based on the decomposition of an arbitrary vector space 
into a direct sum of two subspaces, the two corresponding orthogonal projectors allowing one 
to define an i?-matrix of a universal XX model, and then of a Hubbard model using a Shastry 
type construction. The construction is very general, since the two XX models that are used 
can be based on different (super) algebras. The QISM approach ensures the integrability of 
the models, leading to local Hubbard-like Hamiltonians thanks to the appropriate behavior 
of the obtained /^-matrices. 

Recent studies of the Hubbard model and its generalizations, were motivated by its 
"recent" applications in = 4 super Yang-Mills theory (SYM) - for example, see [18] 
and references therein. Although it appeared that the one-dimensional Hubbard model is 
not the proper answer due to transcendental contributions to the anomalous dimensions (in 
the su{2) subsector of the theory), one may find new directions in this field by studying 
integrable extensions of the one- dimensional Hubbard model. Another aspect lies in the 
possibility of applications to condensed matter physics, particularly when dealing with small 
rank algebras. 

In this paper, we focus on a subcase of the "universal" presentation where the superalge- 
bra for one XX model is gl(n\m) while the other XX model is based on the algebra gl{2), the 
projectors for both XX models being one-dimensional. This kind of model is simple enough 
to investigate further the problem of Bethe Ansatz Equations (see [17]), and presents in- 
teresting features to find some applications in condensed matter physics as well as some 
relevence in AdS/CFT correspondence. 

The paper is organised as follows. In section [2] we prepare the notations and recall some 
results on generalized Hubbard model based on unitary (super) algebras. In section E], we 
solve via the coordinate Bethe ansatz a toy model constructed on gl{2\l) © gl(2) with one- 
dimensional projectors and we generalize the obtained results to gl{n\m) (B gl{2) models in 
section m In section [5l we study a model based on gl{2\2) © gl{2), where the projector for 
the gl(2\2) part is now two-dimensional. We argue on the relevence of this model for the 
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AdS/CFT correspondence. In section [6|, we perform the Jordan- Wigner transformation for 
gl{2\l) © (7/(2) and gl{2\2) © gl{2) models to construct physical Hamiltonians. 



2 Some results on ^/(n|m) Hubbard model 

We remind in this section the necessary notations for self-consistency and understanding of 
the paper. For more details, we refer the reader to [17]. 

The starting point is the definition of the i?-matrix of an XX model based on the super- 
algebra Q = gl{n\m.). Consider the graded vector space V = C"!"" (with possibly m = 0) and 
define the Z2-grading on indices j by 

11 tor n + 1 < J < n + m. 

Accordingly, the elementary matrices E^^ (with entry 1 at row i and column j and else- 
where) have grade [E^^] = [i] + [j]. We introduce a set of integers Af G {l,...,n + m} and 
AT its complementary subset in {1, . . . , n + m}. We then define the following projectors: 

7r = J2 ' ^ = ^^^^ = (2-2) 
as well as (using auxiliary space notation) 

S12 = TTi 7f2 + Vfi 7r2 . (2.3) 

The i?-matrix of a XX model based on (7/(n|m) with projectors {n,W) is given by: 

R12W = S12 P12 + S12 sin A + (I ® I - S12) P12 cos A (2.4) 

where P12 is the graded permutation operator and A G C the spectral parameter. It obeys 
Yang-Baxter equation, is unitary and regular. 

The i?-matrix for (generalized) Hubbard models is obtained by coupling the i?-matrices 
i?|2(A) and Ri2W of two independent XX models, the coupling constant being related to 
the potential U of the Hubbard model under consideration [8,9, 15-17]. We stress that the 
two XX models can be based on two different (super) algebras and Qi, with two different 
(graded) vector spaces V| and and two different projectors and Tr^, associated to two 
different sets A/| and A/j. 

Introducing the parity matrix Ca (« =T or i)- 

C^ = J2 - E = - ' (2.5) 

the i?- matrix of the Hubbard model based on the pair of (super) algebras Q-^ and Qi is given 
by: 

RliiX,, A2) = i?l2(Ai2) R^Xu) + tanh{h[,) RlM,) RiM,) (2.6) 



2 



where A12 = Ai — A2 and X'12 = Ai + A2. In the same way, = h{Xi) + h{X2) and the 
function h{X) is such that 

sinh(2/i) = U sin(2A) . (2.7) 
Note that the site 1 for is composed from the tensor product of the site " 1 | " appearing 

T I 

in the matrix by the site "1 i " which is in the matrix i?{2- This is obviously the same 
for any site we will consider in the following. 

The i?-matrix (12. 6p is symmetric, regular and satisfies the unitary relation. Moreover, 
when the relation (12. 7p holds, the i?-matrix (12. 6p satisfies the Yang-Baxter equation: 

Riii-^i, ^2) RiiiXi, X3) rI^{X2, X3) = rI^{X2, X3) rI^{Xi, x-s) r12{Xi, X2) ■ (2.8) 

Being equipped with an i?-matrix with all required properties, we can proceed to define the 
corresponding quantum integrable system, by performing the following steps: monodromy 
matrix, transfer matrix and Hamiltonian. The L-site monodromy matrix is given 

T,<t....t,>iX) = RlliX,0)...RlliX,0) (2.9) 

and its transfer matrix is the (super)trace in the auxiliary space: 

t(A) =tr,T,<6,,..,,>(A). (2.10) 

Then the generalized Hubbard Hamiltonian reads 



H^^^lntiX) 



= J2Hx,x+i (2.11) 

A = a;=l 



with 

Hx,X + l = (SP) + (SP) + uC^kCik, (2.12) 

where we have used periodic boundary conditions. The notation x,x+i means that the 
operator O acts non trivially in the parts x | and (a; + 1) t only. It acts as identity on all 
the sites different from x and a; + 1 and also on the parts x | and (a; + 1) | of sites x and 
X + 1. Explicitly, one has 

(sp)..,.+i = + 

with a = t or |. The indices j and J run over Af and 77 respectively. 
One can also define the momentum operator: 

exp(2p) = t(0) = P,L P2L--- Pl^i,l ■ (2.13) 



3 gl{2\l) e gl{2) model 
3.1 Preliminaries 

We consider an example of the above model for a particular algebra gl{2\l)^ © gK'^)l- This 
notation means that for different spin value up or down we take different algebras in the 
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construction of the Hamiltonian. The generic expression of the L-site Hamiltonian is still 
the same: 

L L 

Hgl{2\l)<S>gl{2) = ^ [(S-P)t x,x+1 + (SP)| x,x+l] + U ^ {C^x C[x) , (3.1) 
x=l x=l 

but we now choose the projectors n-^ and tt^ to be such that 

/V P\ Z?12 7-i21 I 7-i21 7-1I2 7-1I3 rp31 , ttiSI ttiIS {o o\ 

(SP)| 3;,^.+! = El^Ef^^^ + Ef^Ej^^+i , (3.3) 

n — ??22 7-133 r' — ??22 /o/|\ 

Performing the Jordan-Wigner transformation one can write the Hamiltonian in terms of 
fermionic creation and annihilation operators (for more details see section [6]): 

L L 

Hgi(2\i)(BgH2) = {cix+iCax + cl^Cax+i] +u^{l-2nl^){l-2n''^^) 

^=1 x=l 
L 

+ ^ |(c|a:+l'^Ta: + x'^l x+l) {^'^ x^%+l " x ~ x+l) 
x=l 

+ c| C| X d\ ^j^i d^x + c\^ci^ x+1 d\ ^ d^ x+i | 

L 

-w^(l-2nJJ(l-n^>?, (3.5) 

x=l 

where = bj^^bax is the particle number operator for b = c,d, and we assume standard 
relations between the operators (a,/3 =1, i): 

'^f3y} — ^xy Sap i {d^x^ df^y} — 5xy 5af3 ] {Cax, dpy} — . (3.6) 

In the next subsection we present the system of Bethe equations corresponding to this 
Hamiltonian. Then, in the following subsections, we explain in details the method we used 
for our calculations: in subsection 13. 31 we give the description of the coordinate Bethe ansatz 
approach for the first step of "nested" gl{2\l) ® gl{2) problem. Next subsections consist the 
explanations of the second and third steps of the problem. 



3.2 Result for gl{2\l) gl{2) 

As we shall see in the construction of the Bethe ansatz, the model describes three different 
kinds of "particles" {e^\ e^^ and e^^) above a "vacuum". The particles e^^ and e^^ can be 
associated to spin up and down electrons, while e^'^ represents a " spin 0" fermion (see section 
|6]for more details). 

The energy of a state with excitations is given by 

AT 

E = L-2N + 2^coskm (3.7) 
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and its momentum reads 

N 

P = ^km (3. 

m=l 

where the "impulsions" (or Bethe paramaters) km obey the Bethe equations of our gl{2\l) 
gl{2) modeh 

^ i sin kj + iam + 7 



^ jj ^ j = l,...,iV (3.9) 
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N ■ ■ 1 , ■ , u K ■ ■ , u 

I sm Kj + + f . , ^, T-r ^c^TO — ^a^ + f 



-1)^TT , ^ ^ =A(r?) TT " - ^ 2 ^ m = l,...,i^ (3.10) 



i=l "* 4 1=1, l^m 

2iTV I 



A(n) = e^("l+■■•+""^ M = 0,...,K, l<ni<n2< ... <nM <K (3.11) 

where L is the number of sites considered in Hubbard model, N is total number of e^^, e^^ and 
e'^^ "particles". K counts the total number of excitations e^'^ and e^^^ and finaly M numbers 
e'^^ "particles". The integers Uj correspond to the Bethe parameter of the last level, but their 
Bethe equation has been solved: it just corresponds to the quantization of these parameters 
(see details in section 14.41 that deals with a more general case) . 

3.3 Coordinate Bethe ansatz, level one 

We use the coordinate Bethe ansatz to find the eigenvalues and eigenvectors of this model. In 
this model we have 5 types of distinct excitations noted e^^, e^^ and e^-*-, e^K The objects 
g^j^g elementary vectors (with 1 in a position and elsewhere) which form a natural basis 
for the vector space Vax ~ (a and x hold to describe spin and position in physical space). 
First of all we define the reference state (pseudo- vacuum) : 

L 

h = X{e':e]}. (3.12) 

x=l 

In the following calculations we use the expression of the Hamiltonian only in matrix 
form (13. ip . The corresponding eigenvalue is given by: 

H gl{2\l)®gl{2) 00 = uLcpQ. (3.13) 

We define an excitation by a pair {A, a), A = 2 or 3 (corresponding to vectors or e'^) 
and a =1, | and there is no pair (3, |). The excitation states of the Hamiltonian can be 
written as 

0[(A a)] = J2 ^[x, (AM e^l -^ . . e^- (3.14) 

X 

where the sum over x = (xi,a;2, . . . ,xjv) is considered without coinciding points xi = Xm 
such that ai = am for any 1,771 = 1, . . . , N (this describes the exclusion principle for identical 
particles). We also use the notation {A, a) = (Ai, ai) . . . (A^v, oat), and we omitted the 
'empty sites', i.e. sites carrying vectors e^^e^K 
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Applying the Hamiltonian to fl3.14p . one gets the eigenvalue equation for the function, 
given by 

N 

J2 (^[x - em, (A a)] A- + ^[x + e„, (A,a)]A+) + 

m=l 

+ - 2N) + uJ2^i^i - ^n)S{ai ^ a„) - e)^[x, (Z,a)]A=^ = (3.15) 

l,n 

where is an elementary vector in with entry 1 on the m*'* position and elsewhere. 
Also we denoted 

J] J]5l(x, ^x„) l[6Hxiy^xJ6Hxiy^Xm±l), (3.16) 

Ij^m n^m I 

HdHxi^Xn), (3.17) 

1 - 6{xi - Xn)S{ai - an) . (3.18) 

All these symbols mean that there is no particle with the same spin on the same and neigh- 
bouring sites with some conditions corresponding to each symbol (exclusion principle). 

It is convenient in the following to denote \l/[x, (A, a)] simply by \E'(x) when there is no 
ambiguity. 

At the first step, we take a non-interacting regime, which implies xi -C X2 ^ . . . ^ xat (in 
other words, the particles are far enough from each other). As a consequence, all symbols in 
fl3.18p are equal to 1 and there is no interaction term. We look for a solution of this equation 
in a form of "free particles", namely ^E'(x) oc e***^^, where {ki, k2, ■ ■ ■ , kj^} denote a set of 
unequal numbers. 

We get therefore the value of energy 

N 

E = L-2N + 2^coskm. (3.19) 

■m=l 

In order to consider all other cases of the particle dispositions, we assume the Bethe 
hypothesis for the general solution of \E'(x). We divide the coordinate space {xi,X2, ■ ■ ■ ,xn) 
into A^! sectors: for Xq^^ < Xq^ < . . . < Xq^^, 

^q(x) = ^(-1)[^''(«)]$(P, gp-i)e*^''-'3- (3.20) 
p 

where P = [pi,p2, ■ ■ ■ ,Pn] and Q = [qi,q2, . . . ,qN] are two elements of the permutation 
group ©AT and Pk ■ Qx = '^^kp-Xq.. The symbol (— 1)['*9('3)1 stands for the signature of 
the Q-permutation when restricted to grade 1 particles (that is e'^^). For instance, we have 
the property (valid for any permutation Q and any permutation Iljj+i): (^—I'^^'^diQ^u+i)] = 
We recall that \E'(x) and accordingly $(P, QP~^) depend on the type of 

excitations [A, a). 

The coefficients $(P, QP~^) are not all independent and the application of the Hamilto- 
nian represented in equation (13.151) can reduce their number in several cases. 



m 



A3 

5^ixi 7^ Xn) 



6 



1. Let us consider the sector Xq-^ <^ . . . <^ Xg^ < a^g^^^ <^ Xg^^ with Xg. = x^.^-^ — 1 and 
ttg. = aqi+i- In this case A~^^ = and A+ = 0, all other symbols in fl3.18p are equal to 1. 
Thus, equation (13.151) becomes 



J2 (^Q(x-e„,) + ^Q(x + e™)) +^Q(x-egJ + ^Q(x + eg^^J = {E - L + 2N)^q{^) 

(3.21) 



Using (13.201) and performing simple algebraic calculation we obtain 

$(n,feP, QP~'Uab) = -$(P, QP~'), (3.22) 



where Uab is the permutation of objects a and b which are linked to i,i + 1 by p~^{a) = i 
and p~^{b) = i + 1. This relation holds for any value of Ag. and 

2. Now let us consider another case: -C . . . <C = x^.^-^ -C Xg^ and ag. ^ ag,+i- We 
denote by Q the sector where i and i + 1 are permuted (i.e. Q = QUu+i). It implies that 
the equation (13.151) becomes 

J2 [*q(x - e„,) + ^'q(x + e^)] + *q(x - BgJ + ^'q(x + e^J + 

+ *q(x - eg^^J + v&q(x + eg^^ J + ((L - 2N) - E) v1>q(x) = . (3.23) 
On the other hand, we have the condition of continuity of the wave function 

*q(x) = *q(x) with Xg, =Xg,^, (3.24) 

These last two equations together give two conditions on the coefficients ^{P,QP^^) which 
can be written in matrix form 

^{YiabP^QP-^Iiab)) \rab U) \^{P, Q P-^Tiah) ) ^ ' 

with 

_ 2i{\a - Ab) _ -u 

u + 2i{K - h) ' u + 2i{K - h) ■ ^ ' 

These equations also hold for any type of excitations (any value of Ag^ and Ag^^J. We set 
\a = sin ka to simplify the expressions. 

Before applying the periodic conditions on the coefficients $(P, Q-P"^) (henceforth we 
will rename QP~^ simply by Q), we rewrite equations (13.221) and (13.251) in a compact form. 
Let us introduce the vector 

|,(P)= <^[P,Q,{A,a)]\{Ag„agJ,...,{Ag^,ag^)> (3.27) 

Q,(A,a) 

where the sum is over all types of excitations and all corresponding sectors. The vector 
|(ylg^, agj, . . . , (Ag^, Ogj^) > belougs to Vi (g) . . . ® Vn, where V = span{2 t, 2 i, 3 t}, and 
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represents one state in the space of excitations. Let us illustrate by an example for N = 2 
excitations: in this case we have 6 types of different excitations. 

J2 ^P^Q^ (A^.aJ >= HP,^d, 1)|2 T,2 T> +HP,td,2)\2 t,2 i> 

Q,(A.a) 

+$(P,2c/,3)|2 t,3 T> +<f(P,ni2,2)|2 i,2 t> +^P,td,A)\2 1,2 [> + 
+$(P,2t/,5)|2 i,3 T> +$(P,ni2,3)|3 T,2 t> +$(P,ni2,5)|3 T,2 i> 
+$(P,zrf,6)|3 T,3 T> • (3.28) 

Thus, for N = 2 excitations we can express fl3.22p and (13.251) by 

$(ni2P) = (-1)$(P, td, 1)|2 T, 2 T> + [ti2<f (P, td, 2) + ri2$(P, Ui2, 2)] |2 t, 2 i> 
+ [ti2$(P, ni2, 2) + ri2$(P, td, 2)] |2 i, 2 T> +(-l)$(P, td, 3) |2 t, 3 T> 
+ (-l)$(P,ni2,3)|3 T,2 T> +. . . + (-l)$(P,zrf,6)|3 T,3 T> 
= S[l\X^-\2mP) (3.29) 



where we introduced the so-called S'-matrix S'jg (A): 





/|2t,2T>\ 
|2T,2i> 
|2 T,3 T> 


/ -1 




ri2 


\ 

-1 


/|2t,2T>\ 
|2T,2i> 

|2 T,3 T> 


*S'l2'*(Al — A2) 


|2i,2T> 
|2i,2i> = 

|2i,3T> 


ri2 


tl2 


-1 

tl2 


ri2 


|2i,2T> 
|2i,2i> 

|2i,3T> 




|3 t,2 T> 
|3T,2i> 

V|3T,3T>/ 


-1 

v 


ri2 




tl2 

-1 / 


|3 T,2 T> 

|3T,2i> 

Vl3T,3T>/ 



(3.30) 



The coefficients ti2 and ri2 are defined in (13.260 . The same reasoning could be repeated for 
an arbitrary number of excitations, and we obtain 

<l>(n„,P) = 5i;)(A, - A,)<I>(P), (3.31) 

where the matrix >S'^),'*(Aa — A;,) acts nontrivialy only on vector spaces Va ® Vf,. It satisfies the 
Yang-Baxter equation: 

'5'l2'*('^l ~ A2)5'|3''(Ai — A3)S'23''(A2 — A3) = 5'23''(A2 — A3)S'j3''(Ai — A3)S'|2'*(Ai — A2) . (3.32) 

3. In order to obtain the Bethe equations we should apply the periodic boundary condi- 
tions on the function ^E'(x). Let C be the cyclic permutation given by C = IIjvi . . . IlTVAr-i. 
More precisely the periodicity condition means 

vI/QC7(x + e,,L) = M/Q(x) (3.33) 

which implies a condition on the coefficients $(P, Q), namely 

|)(PC) = e*'=w^<|.(p) (3.34) 
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choosing P = ^ , one derives a system of equations satisfied by the coefficients ^{id) 
which is called "auxiliary problem". For j = 1, . . . , A^, it reads 

S^%^ . . . S^i]s[f . . . Sf\^^td) = e'^^'^^td) (3.35) 

here we omitted the arguments in S'-matrices, S'^f^ = S''^fj{Xa — Xb)- 

If we could perform the diagonalization of the left-hand-side in the general case, it would 
be possible right now to write the Bethe equations for our model. The form of matrix S^,^ 
requires to use again the coordinate Bethe ansatz to solve the auxiliary problem. 



3.4 Auxiliary problem, level two 

Again we have the equation on eigenvectors and eigenvalues to solve, but contrarily to the 
first step the Hamiltonian is a little more complicated: 

S!^,^...S^S^...S^,^<P = A,<P (3.36) 

Since the matrix S^^^ is regular, the above Hamiltonian can be identified with t{Xj), where 
we have introduced a new transfer matrix 

t{X) = tro ^5']^\ q(Aj+i - A) ... S^j^^Q^Xjy - A) SjloHAi - A) ... 5]^^ o(Aj_i - A) o (Aj - A) j 

Since S^^^ obeys Yang-Baxter equation, this new problem is indeed integrable. We can thus 
perform the Bethe ansatz again. 

First, we slightly modify the S'-matrix = 5*12 to simplify the following 

calculations. Thus, it is given by 

^(Ai - A2 
^(Ai-A2) + f 

= J2 (3-38) 

i,i=l,2,3 
1=2,3 

We use the coordinate Bethe ansatz to proceed with this problem. At this step we have only 
three types of different excitations: e^\e^^ and e'^^ At first, we choose the reference state 
(pseudovacuum) to be 



^i2(Ai-A2) = Pi2-ti2i^i2, With = TTT^^xT^ ' (3-37) 



TV 



k=l 



with eigenvalue A = 1. 

Applying the S'-matrix on different states (on the pseudovacuum and two excitations 
and e^), we get 

^i2(Ai - A2) ® e^i = e^i ® e^^ (3.40) 

5i2(Ai - A2) ® e'^ = -ti2 ® e'^ + (1 - h^) e'^ ® e^, (3.41) 

Si2(Ai - A2) ® = -tu e'^ ® + (1 - ^12) ® e^^ (3.42) 

5i2(Ai - A2) ® = (g) with A, B = 2 ^ 3 ] . (3.43) 
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We first consider tlie "non-interacting" regime, which imphes to solve the equation with only 
one excitation. Let A = 2 f or 3 t and indicate the type of the excitation 

0(A) = ^ f^^A . (3.44) 

X 

We will forget henceforth the A index to simplify the notations. 

Here, contrarily to the first subsection, it is not obvious to readily calculate the Hamilto- 
nian action on excited states. We can do it step by step introducing a recursive Hamiltonian 
(for details, see [19]) 

Sj^kj ■ ■ ■ Sj-ijcj) = (3.45) 

(k) 

and recursive coefficients fx which represent the action of the recursive Hamiltonian on the 
coefficients fx- One gets the following relations 



: (3.46) 

fi'^^'^x = ^j-kj ■ ■ ■ Sj-lj ^ fxG^^ ■ 



We write down the recursive relations between coefficients /i'^'* and /i*^ ^\ For x j,j — k 

jik) ^ j(k-l) ^ 

= -V../r'^ + (l-V..)/£^ (3.47) 



Ij ''3-kjJj ~ V-^ f-j-kjjjj-k 

dk) . . -f C^-i) _L n _ + . ..^fC"-^) 



fj-k - -h-kjfj-k + (1 - 'tj-kj)fj 

Skipping the details of calculation, the solution of the equation (13.451) gives a relation between 
different f^. Closer considering shows that one needs to introduce an additional constant a 
(particle rapidity), such that 



fx+i{a) iXx + ia + ^ 



(3.48) 



fxia) iXx+i+ia-^ 

Thus, we can write down the expression for coefficient fx and the eigenvalue corresponding 
to one excitation. We normalize /i = 1, then we have 

m=l ^ ' 

A- = a J a) = — f . 3.50 

(k) 

It is convenient to write down some relation between iterated coefficient fx , fx and a: 

fx{a) = ai{a) . ..ax-i{a) 



iXi + ia - f 

iXx + ia — ^ ^ 
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/f(a)= //""^ (3.51) 
fj%{a) =aj{a)fj_k{a). 
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Starting from the initial equation (13.361) . we can choose j = N and without any loss of 
generality we can derive the equation which represents the Bethe equations on rapidity a by 
application of the initial Hamiltonian on /^(a). It also automatically implies the periodic 
boundary condition 

N 

Y[am{a) = l. (3.52) 

m=l 

All these results hold for any type of excitation e^^ or e^^ (we neglected the index A above). 
Now we can proceed to the case of K excitations. Let A be a vector {Ai, . . . , Ak) with 
= 2 I, 3 I, which represents the combination of different excitations in consideration. The 

eigenvector is naturally represented as 

^{A) = J2'^{^,A)e^^...e^-. (3.53) 

X 

The sum is done over all coordinates Xi without coinciding points Xi = Xk for any i, k 
(exclusion principle). 

We use the Bethe hypothesis for the coefficients \E'(x, A), in the sector Xg-^ < Xg^ < . . . < 
Xqj^ where Q = [gi, . . . , g^] is a permutation of the integers 1, 2, . . . , /T: 

^(x, A) = Y, ^{P. QP-\ ^)hq, MU, M ■ ■ ■ (a,,) (3.54) 
p 

where /x, (opj is the one-particle solution with rapidity Op.. We will omit the index A to 
simplify the notations. As above, one can reduce the number of coefficients $(P, QP~^). 

At this stage we will consider K excitations of different types. Therefore we have 
K\/M\{K — M)\ sectors x^^ < Xg^ < . . . < x^^, where M, K — M are the numbers of 
excitation of type 2 | and 3 | respectively. Acting with the recursive Hamiltonian on the 
eigenvector 0, see (13. 45 p . we are able to write down the relations between different \E'q''(x). 

We consider j to take the biggest value such that we do not have any excitation on the 
sites bigger than j. Several cases can occur, depending on whether there is an excitation on 
sites j and/or j — k. If there is an excitation on site j, it must correspond to Xg^ = j. For 
the site j — k, there exists such that we have either Xg^ < j — k < (no excitation 

at site j — k), or Xg„^ = j — k (one excitation at site j — k). Then, we have 

vl>(')(x) = vl>g-^)(x), for X ^ j,j - k, (3.55) 



■^QK IK IK 

^^q\- ••,],•••) = -t,~-k,^V\- • • , j , • • .) + (1 - t,^k,)^^^r'\. . . , J - fc, . . .), 

with Q' = QIlKra+i ■ ■ ■ for x=^j-k, 



-^Qm -^Qm •'^qm 

with Q' = QUmK ■ ■ ■ IImm+1 

for X 7^ j 

In this last case, we have also the exclusion principle: for any n ^ m, Xn ^ j,j — k. 
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When two excitations are on sites j and j — k, with = j and Xq^ = j — k, we have: 

^ (. . . , J — fc, . . . , J, . . .) = ^Qn^^^ {■■■■, j-, j — k, .. .) for non-identical excitations, 
(. . . , j — A;, . . . , j, . . .) = \I/q ^^(. . . , j — /c, . . . , j, . . .) if the excitations are identical 

On the one hand, we have the Bethe hypothesis for the coefficients ^'(x) composed of the one- 
excitation functions fx{i on the other hand, we have the relations for the iterated coefficients, 
it is natural to calculate how the iteration passes to the one-excitation functions. In other 
words, we want to calculate the coefficients '^^q (x) in terms of the iterated free excitation 

functions f^\ap). 

After first iteration for any x except the coefficient • • , J — 1, • • • • • •) we find 

W = E *(^' QP^")f^i Mfill M ■ ■ ■ ffl M ■ (3.56) 

p 

The next iteration will allows us to make an assumption on the form of \E'q''(. — 

1, . . . , j, . . .). Thus we calculate the coefficients \1/q''(x) except again those where two co- 
ordinates coincide with j — 2 and j. Skiping the calculations we find for the first iteration 

■ ■ , J - 1, • • • , J, • • •) = mKK-i)P-')f^lM . . . ff\a,,^:)ff\{a,,). 

p 

This implies the conditions 

^q'*(. . . , j — 1, j, . . .) = \EfQ(. . . , j — 1, j, . . .) for identical excitations (3.57) 
. . , j — 1, j, . . .) = '^Q'{- . . , j, j — 1, . . .) for different type of excitations (3.58) 

with Q' = QUkk-i- It is useful to notice that for identical excitations QUkk-i defines the 
same sector as Q. These equations give the connection between ^{PUkk-i, Q{PIIkk-i)~^) 
and ^{P,QP"^) for some sectors Q. We will consider it precisely in the general case. 
For the second iteration we find 

W = E '^(^' M ■ ■ ■ fSl M (3-59) 

p 

except 

• • , J - 1, • • • , J, • • •) = E '^(^' mKK-i)P^')ffl{a,:) . . . ff\a,,_,)ff\{a,,). 

p 

(3.60) 

and as in the previous case we have two "undefined" coefficients ^!^'^\. . . , j — 2, . . . , j, . . .) 
and ^(2)(...,j-2,...,j-l,...,j,...). 

Generalizing these results until the {k+l)-th iteration, we make again a set of assumptions 
for the coefficients of type \l/*^'^^(. . . , j — k, . . . j, . . .). Hence, for the k-th. iteration, we find 
the coefficients 

^?(x) = E •^(^' • • • fit M (3.61) 

p 
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except a set of coefficients. Let rii be integers such that < ni < n2 < ■ ■ ■ < rii < k, 

I = 1, . . . , k — 1, that label the position of the possible excitations between j — k and j. We 
have a set of assumptions made from (k + l)-th iteratiorH. 
If there is an excitation in position j — k, we have: 

^qU - = E •^(^' (Qnf-i)p-i)/W(a,J . . . /f (apK_J/j-U«.J (3-62) 
p 

m^^\j-k,j-n,,j)= (3.63) 



■^^Q\3-k,3-nu...,3-n^,3)= (3.64) 
m (Qnf-'-i . . . nf-^)p-i)/W(a,J . . . /f (apx-,-J/j-U«PK-J • • • fflMvK) 



^S^(j-A;,j-/c + l...,j-l,j)= (3.65) 

Y ^{p. mK-' ■ ■ ■ nf-^)p-^)/i:j(apj . . . /f Vp.-j/j-U«PK-..j • • • /f 

p 

If there is no x equal to j — k, from the calculation we have 

- ^1, j) = E *(^' (^?nf-^)P-^)/itJ(a,J . . . /f VpK-J/f (3.66) 
p 

*Q^(i -^^2,7 -ni,j) = 

Y HP, mr'^f<~')P~')fi':lM . . . /f (SK-J/A,(ap.-J/i-l,(sJ, (3.67) 
p 

vi>(J)(j-A: + l,...,j-l,j)= (3.68) 

X: $(P, (Qnf-'=+^ . . . nf-i)p-i)/W(a,J . . . /f (a..-.,J/j-Ui(«PK-..J • • • /i-U«.J 
p 

We can write a set of conditions from all made assumptions to determine the relations 
between different $(P, P~^). One can see that we can consider only the fc-th iteration to 
derive the necessary relations. More exactly we have k equations to satisfy 

^Q^(i - k,j -ni,.. . , j - ni,j) = '^Q~^\j - k,j -ni,.. .,j - nij) (3.69) 

for identical excitations corresponding to j and j — k. Remark that in this case QH^"'"^ is 
identical to Q. 



^In order to make the notation shorter we do not write the dots . . . for unchanged indices when it is not 
ambiguous. 
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For non-identical excitations, we have 

^Q^(j -kj -ni,...J - nij) = "^Q^K-i-iUJ -ni,...J- nij - k) (3.70) 
with / = 0,...,/c — lin both cases. Hence we have a set of conditions: 
Ep hp, (Qnf-'-i . . . nf-i)P"i)/w . . . f^'\a,,_,_Jfl%{a,,J . . . (a,J = 

= Ep m m^-'-' ■ ■ ■ nf -^)p-i)/('=-i) . . . j/f -'Vp.-J • • • ft^M ■ 

Using the expressions obtained for iterated one-particle functions (I3.5ip and omitting all 
the algebraic calculations, we have the expression valid for any sector Q and any type of 
excitations: 

$(pn,+i,Q(pn,+0-^) ^ "^"^^^ = -^sL - f ^ = ^-^,---,^-i- (3.71) 

Before applying the whole Hamiltonian (13.361) . which also represents the periodic bound- 
ary conditions, on the coefficients $(P, QP^^), we rewrite the obtained equation (13.711) in 
a form gathering all possibilities of excitations A = [Ai, . . . , A^). At first we introduce an 
object (henceforth we will rename QP~^ simply by Q) 

$(P) ^ ^ $(p, :4) 1^^^, ...,A,^> (3.72) 

where the sum is over all types of excitations and all corresponding sectors. The vector 
\Aq-^, . . . , Agj^ > belongs to Vi ® . . . ® (where V = spanje^^, e^^}) and represents one 
combination of K excitations. To explain the notation we write an example for K = 2 
excitations. In this case we have 4 types of different states: 

ZqaHP,Q,^)\K^K >= HP,td,l)\2 t,2 T> +<f(P,2rf,2)|2 T,3 T> + 

$(P,ni2,2)|3 T,2 T> +<^>(P,^rf,3)|3 T,3 T> • 

Thus, for K = 2 excitations we can express (I3.7ip by 

HUuP) = a^^' (^HP.td, 1)\2 t,2 t> +<f(P,2rf,2)|3 T,2 T> + 

$(p,ni2,2)|2 r,3 T> +$(P,^rf,3)|3 T,3 T> ) = ^!?<^'(^) 



(2) 

where we introduced, as in the first section, a new "^-matrix", noted 5*^2 • One can see that 
this S'-matrix in this case is a pure permutation. 

The same reasoning could be repeated for an arbitrary number K of excitations and we 
obtain 

mabP) = C^JPaMP), (3.73) 

where a = p{i), b = p{i + 1) and the permutation Pab acts nontrivialy only on ® H vector 

(2) 

spaces. Obviously, the matrix being a permutation, satisfies the Yang-Baxter equation. 
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To obtain the Bethe equations we apply the periodic boundary conditions which are 
automatically implied if we pass to the initial equation (13.361) . We choose j = N, to get 

If no X equals A^, we find the eigenvalue An = crAr(ai) . . . o"Ar(ax) where a is the eigenvalue 
of one excitation function (13.501) . Otherwise, the equation 

^i^~^\xu ...,N,...,xk) = An^q{xu ...,N,...,xk) 
leads to additional condition on the coefficients Q). Hence, in terms of ^{P), we have 

N 



k=l 



where C is a cyclic permutation given by C = n^i • • • ^kk-i- 

Choosing P = C^~™, we can derive a system of equations on coefficients ^{id) which 
will be called " auxiliary problem 2" : 



N 

5iV---4l5S---5m-im^(^^) = n^'^(«-)'^(^^) m = l,...,K. (3.74) 



k=l 



At this stage, we arrived to the third step of nested coordinate Bethe ansatz. We have 
only 2 types of possible excitations e^^ and e'^^ and the S'-matrix (I3.73P is a pure permutation. 

3.5 Permutation problem, level three 

As mentioned above, we have only 2 types of "particles": e"^^ and e^^^ and an equation 
composed only with permutations: 

n, 4> = Pj+^,...PKjPij...P,-ij <i) = A(t). (3.75) 

Again, we solve the problem with the coordinate Bethe ansatz, and fix the first "particle" 
e^^ as the "new" vacuum state 



K 

>M=Q 



Hef, A = l. (3.76) 



i=l 

The second "particle" (e^^) represents some excitations above this vacuum state. One can 
verify that the function 

M 

0M(n)= ^(^)n^S (3.77) 

with 

M 

v|/(f) = '^iP)X{9.A(^P(^)) , Ua) = (3.78) 

P&Sm *=1 

15 



is an eigenfunction of this permutation problem if the coefficients $(-P) obey some conditions 
that we formulate below. 

Acting with the Hamiltonian Tij gives only the periodicity condition on the coefficients 
$(P): 

^(^p^ ^ I'^pw]^ ' ^ ^ Him-.-Hm-im • (3.79) 
We assume a simple condition on the coefficients $(P): 

$(Pn,,+i) = $(P) for ^ = 1,...,M. (3.80) 
Imposing these conditions and periodicity leads to 

(7.(a(n)) = e'F"^ n = l,...,K. (3.81) 
Therefore the eigenvector can be written in the form 

M 

(f)M{n) = Yl XI n^^"''''''^'^S wi^h 1 <ni <n2 < ... <nM < K . (3.82) 

xi<X2<...<xm P&Sm *=1 

It has the eigenvalue 

A(n) = fje"^"^ = exp ('^ \n\] . (3.83) 
j=i \ / 

Then, the total number of states (pM{ni, ....um) is 

5: ' - ( M ) = mkk'^ My. 

which is the right number of eigenfunctions with M excitations. Hence the conditions (13.801) 
and the functions (13.811) solve the permutation problem (13.751) . 

Gathering the results of sections 13. 3i 13.41 and 13. 5^ we get the Bethe equations written in 
section I 



4 Generalization to ^/(n|m) 5'/(2) model 

In this section we generalize the results obtained in previous sections to larger algebras, 
namely the case of (yf/(n|m)| © 5'^(2)j Hubbard model. The Hamiltonian is given by the 
expression 

L L 

Y [i^Ph + (SP)i x,x+i] + xC^ x) , (4.1) 

x=l x=l 

with the notation: 

ti+m 

(SP)t.,.h-i = E(^Tx^?".+i + (-l)^'^'^|x<.+i) (4-2) 

a=2 



imix,x-,i = EllElU^ + EflEll^, (4.3) 

n+m 

C^x = - Y ' ^l^ = Eli - Ef^ (4.4) 



a=2 
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and the grading we use is given in (12. ip . 



4.1 Result for gl{n\m) gl{2) 

We first gather all results detailed in the following sections and write down the Bethe equa- 
tions of our model gl{n\ra) © gl{2): 

^(_^)X+A.+i-Q ^sinfc,+za^ + t ^ j = l,...,7V (4.5) 
^ I sm kj + lam - I 

-|-^ ^sinfc,+2a^ + ^ ^ „ ^"--^"' + 2 ^ = i,...,ir (4.6) 

M 

?(3)N / f^|-7(3)| \ 1^(3), V-^(3) 



A(n(3)) = exp — , = 



0<M <K and 1 < nf^ < < ... < 4f < (4.7) 

where L is the number of sites considered in Hubbard model, N is the total number of e'^^,e'^^ , 
e^^,...,e*^""'"™)^ "particles". K counts the total number of excitations from e"^^ to e^"'^'^^^ and 
finaly M numbers the e^\...,e^"^"^^^ "particles". 

There are Bethe parameters n[''\ 3 < < m + n, for each particle but they don't 
show up in the Bethe equations. In section 14. 4[ it is shown more precisely how all these 
remaining parameters (that are quantized) appear in the Bethe ansatz construction. 

The energies associated to these states are given by 

AT 



E = {L-2N) +2Y^ cos k„r (4.8) 



m=l 

and their momentum reads 



N 

. (4.9) 

771=1 



Let us note that the Bethe equations for gl{n\m)®gl{2) are very close to the ones obtained 
for 5^/(211) © gl{2). This is due to the particular projectors we have chosen, see eq. (14. 4p . 
More general models (and Bethe equations) can be obtained varying these projectors: we 
come back on this point in section [51 



4.2 Coordinate Bethe ansatz, level 1 

We solve this model via the coordinate Bethe ansatz. In this model we have n + m + 2 types 
of different "particles" denoted by e^^, e^^, e*^""*""^^^ and e^^^e^^. As in previous case we 
choose the vacuum as 

L 
x=l 
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The excitations above the vacuum state are given by 

0[Z] = 5^v[;[x,Z]ef;...ef; (4.10) 

X 

with indices Aj = 2 |, 3 |, (n + m) 1,2 | corresponding to vectors e^^, e'^^,...,e^""^"'^^, e^^. 
The sum over x is again considered without points where two "particles" with the same spin 
are on the same site. describes the number of all excitations and goes from 1 to L. In 
fl4.10p . the sites carrying vectors e^e^K iiot associated to any excitation, have been omitted. 

Now we assume the Bethe ansatz for \l/(x) and follow the steps detailed in section [31 We 
divide the coordinate space (xi, X2, .., xat) into A^! sectors. For Xg-^ < Xg^ < .. < Xg^, we have 

^q(x) = ^(-1)["^'(«)]<I>(P, gp-i)e*^'^-'3" (4.11) 
p 

where P = [pi,P2, --.Pn] and Q = [qi, q2, .., qw] are two permutations of the integers 1, 2, .., 
and Pk ■ Qx = J^i^Pi^qi- -^^ previous section, the symbol (— 1)[*9(<3)] stands for the 
signature of the Q-permutation when restricted to fermionic particles e^'^'^^^^ e*^""'"'"^^. We 
recall that \&(x) and accordingly $(P, QP~^) both depend on the type of excitations A. 
We gather all the coefficients QP^^, A] in a vector 

l>(P) = j2'^[P,QP-\A]\Ag,,..,Ag, > 

Q,A 

where the sum is done over all possible types A and all corresponding sectors Q. A vector 
...,Ag^ > represents one state of A^ possible excitations and belongs to Vi ® ... ® Vn 
(where V = span{2 j; 2 T; 3 j; 4 T; (n + m) T})- 

In the case of only two particles (A^ = 2) we are able to introduce the S-matrix: 

l>(ni2P) = ^;2^(Ai-A2)<l'(P). 

Here, 5'j2'*(Ai — A2) = acts on elementary vectors |yli,yl2 > as 

sS\2i,2i> = -|2i,2i>, (4.12) 
sS\2i,A> = ti2|2i,A>+ri2|A,2i>, (4.13) 
S[^^\A,B> = -\B,A>, for A, P = 2 t, (n + m) t (4.14) 



where 



2i(Ai - A2) -u 

tl2 = — TTTT TT' ^12 - 



n + 2z(Ai-A2) n + 2z(Ai-A2) 

For an arbitrary number of excitation A^ we have 

l>(n„,p) = 5«(A.-A,)$(P), 

where the matrix S'^^^ acts nontrivialy only on Va ® H vector spaces. 
Again, the periodic boundary condition on the function \l/(x) 

^'qC7(x + Bg^L) = ^fQ(x) with C = n^vi-n^VTV-l 

leads to the second step in our problem: 

S^%...S^ji]s[f...S^%^td) = e^'=^^<l(zd) J = 1, . . . , AT . (4.15) 
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4.3 Auxiliary problem, level 2 

As in the previous section, we transform slightly the 5- matrix 5*1^2^ —>■ —Sl^ = S12 to simplify 
the calculations. The idea is that we work in different steps to diagonalize these matrices 
5*12 via the coordinate Bethe ansatz. At each step we specify a "new" vacuum and "new" 
excitations. At this step we have n + m types of different excitations: e^^, e^^, e^^, . . . , e^""*"™^^. 
Choosing the "new" vacuum as 

TV 
k=l 

the state of K excitations of any type is written as 

0(A) = 5^v[;(x,A)e^^..e^- (4.17) 

X 

where the sum is done over all coordinates Xj without coinciding points Xj = Xk for any i, k 
(exclusion principle). A is a vector {Ai, .., Ak) with = 2 | , 3 (n+nr) t (corresponding 
to e^^,e^^ ...e^"^'"''^). Again, sites carrying e^-^ (no excitation) have been omitted in (14.171) . 

The Bethe ansatz for the coefficients \E'(x, A), in the sector Xq^ < Xq^ < ... < Xqj^ where 
Q = [gi, qx] is the permutation of the integers 1, 2, i^, is given by 

v^q(x,A) = Y,^{P,QP-\A)U^^{a,,)U^^{a,,)...U^^{a,,). (4.18) 
p 

fx{(i) is the one-particle solution (for any type of excitation A) with rapidity a: 

fM = n (- .f'^°+i ) . (4.19) 

The eigenvalue corresponding to this state <t>{A) takes the form 

= a,(ai)...a,(ax) (4.20) 
where aj{a) is the eigenvalue of the one-particle solution 

, , i\j + ia + , , 

Next we gather all the coefficients $(P, QP~^, A) in a vector 

$(P) = Yl '^(^' QP'\^M,i^ > (4.22) 

Q,A 

where the sum is over all types of excitations and all corresponding sectors. The vector 
\Aq-^, .., Agj^ > belongs to Vi (8> ... (8> Vk (where V = span{e^^; e^^; e'^"^'"^^}) and represents 
one combination of K excitations. 

Application of the Hamiltonian fl4.15p on the excited state (piA) gives two types of condi- 
tions imposed on the coefficients $(-P). From the first one, in the case of only two "particles" 
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[K = 2), we can introduce again a S-niatrix corresponding to the second step. The second 
type imphes the periodicity condition. Thus, the S'-matrix is defined by 

l>(ni2P) = ^(?(ai - a2)$(P) (4.23) 

with 

(2)/ '^'^1 ~ '^'^2 + f 

Sl^[ai - a2) = ai2-Pi2, and = -. : ^ . (4.24) 

tai — ta2 — 2 

For an arbitrary number K of excitations we have 

^n,,P) = S^\a,-a,MP), (4.25) 

where p~^{j) — p^'^{i) = 1 and the permutation S'^\ai — a-,) acts nontrivialy only onVi®Vj 
vector spaces. The matrix {cii — 0-2) satisfies the Yang-Baxter equation since it is a 
permutation. 

The periodic boundary conditions on $(P) imphed by the action of the Hamiltonian 
f l4.15p is written in following form: 

N 

^Slin.-41^S-^i'ii™<^'M) = n^'^(«-)'^(^^) m = l,...,K, (4.26) 

fc=i 

where the S-matrix arguments were omitted for simplicity. 
4.4 Permutation problem, level 3 

Thus, we arrive to the third step of nested coordinate Bethe ansatz. Here, we have e^^, e^^, ... 
, e''"^"^''^ "particles" that move 'freely', the Hamiltonian T being constructed on permutations 
only: 

T<f) = Pj+i,...PKjPij...Pj-ij = A . (4.27) 

Note that F is a cyclic permutation, and is independent from j. 
We choose the "particle" e^^ as the vacuum state: 

K 

<pM=o = l[ef, A = l, (4.28) 

i=l 

and introduce the function 

M 

<pf^{A) = J2^{x)l[ei:, A, = 3T,4T,...,(n + m)r . (4.29) 

X i=l 

It describes a state with M excitations above the vacuum state 0a/=o- 

The coefficients "^{x) are defined in the sector 1 < a:g(i) < Xg(2) < ... < < K, with 

Q e 6m, by 

M 

^q(^)= E •^^'H^,Q^~')n^-,w(4(l))' 9x{a) = a\ (4.30) 
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/ON 

One can verify that (j)\l{A) is an eigenfunction with the following eigenvalue 



M 



A = n«f (4.31) 



i=l 



if some conditions, which we precise below, are satisfied. 

Application of the Hamiltonian gives only the periodicity condition on the coefficients 

$(3)(p,gp-i): 



$(^)(pc,gp-^) _ (3) 
$(3)(p,gp-i) ^ 



aJiv/-)]^ 5 C — I\.iM---^M-iM ■ (4.32) 



As in section l375l we assume some relations on the coefficients $^'^-'(P, gP~^), but, now, the 
form of these relations depend on whether the particles are identical or not. If, for a given 
i, Xi and Xj+i correspond to identical particles we impose 

$(3)(np(,)p(,+i)P,gp-i) = <|.(3)(p,gp-i), (4.33) 

while, otherwise, we set 

$(3)(np(,)p(,+i)P, gp-inp(,)p(,+i)) = ^^''\p, qp^') . (4.34) 

As we can see, there is a sector changing in the relations above, and we proceed recursively 
using the same methods as above. We introduce 

¥'\P)^Y.^^'\P,Q,A]\A,,,..,A^^ > (4.35) 

Q,A 

where the sum is over possible types A and all corresponding sectors Q G &m. The vector 
\Aq-^, Ag^ > represents one state with M excitations and belongs to Vi ® ... ® Vat (where 
V = spanje'^^, e^^, e*^""'"'"^^}). Then, relations (14.331) and (I4.34p can be rewritten in the 
following form 

l>(3)(n,,P) = Si^^^iP) , = Pat . (4.36) 

Therefore, the periodic boundary conditions on $(P) implied by the action of the chain 
of permutations is written as 

Pm+lm...PMmPlm...Pm-lm^^'\td) = [a^^^ ^^'H^d) , m= 1,...,M (4.37) 

but here we have already only e^^, ... ^e^""'"'"^^ "particles" involved in the calculations. Thus, 
we arrive to the next level of nested coordinate Bethe ansatz, with, again, an Hamiltonian 
built on permutations only, and a new chain of length M. 

Using the previous considerations, we repeat the same method and we "eliminate" one 
by one the "particles" e^^, ... up to e^""'"'""^-*^, choosing it as the vacuum state at each nested 
level. 

We suppose that we have M3 "particles" of type e'^^, M4 of type e^^,..., Mn+m of type 
g(n+m)T^ so that M3 + M4 + ... + M„+TO = M. At each level k = 3, ...,n + m - 1, we have 
particles e'-'^"''^-*^, e*^""'"'"^^ as different types of excitations above the vacuum state built 
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on e'^^. The eigenvector can be written in the same form as in fl4.29p and fl4.30p . 

with the set of Bethe roots {Q('=+i)}A4+i+...+Afn+m ^-^^ coefficients ^''''+^\P, Q p-^) with 
P,Qe &Mk+i+...+M„+„,- These coefficients are used to write the vector ^^''^^\id) that obeys 
the periodicity condition: 

P™+i,„....PM.,,+...+M.,.,™Pi,n...P^-i,„.<l'('=+'n^^^) = [a(^^+i)]^'=+-+^"+"' l>('=+i)(zrf), (4.38) 

for m = 1, (Mfc+i + ... + Mn+J. 

Also, we find the periodicity condition of the previous level (when we pass from level 
A; — 1 to level k) 

Mfc+i+...+M„+„ 
i=l 

form = l,...,(Mfc+i + ... + M„+m) and /c = 3, n + m - 1 with Ms = K - {Ms + ... + Mn+m)- 
At last level, we have only one type of excitations e*^""'"'"^^ on the vacuum state e'-"^'"^^-'^ 
and we can see that it is the same case as in the permutation problem of the model gl{2\l)(B 
gl{2). Thus, using the relation 04.391) for k = n + m — 1, we find the following Bethe equations 
which link the case of one type excitation e^""'"'"-*^ and the previous level with two types of 
excitations {e{"+'"-i)T, e("+'")T].: 



Mn4 

H z = 1, . . . , M„+^_i + A4+^ , (4.40) 

5=1 



together with the result obtained from gl{2\l) © gl{2) model 

(n + m) 

n . 

(n+m) _ 27ri 



^in+-;^e^'"M„+_,+M,+„^ j = l^,,,^M,^m (4.41) 



1 < nS"+'") < n^r""^ < ... < < M„^.^_i + M„+^ . (4.42) 

In the same way, we can write the Bethe equations corresponding to the transition be- 
tween level with ^e^"+^-^)^ ^ g(n+m)T j ^j^^j ^j^g previous one 

A^n + m - 1 + A/n + m 

j^(n+m-2)jM„,„_3+...+A.,„„ ^ -Q ^ = I, . . . , M„+„_2 + ... + M„+^ (4.43) 

and we can continue this recurrence up to a^^\ 

Solutions for every set of Bethe roots a'-'^^ can be computed as following. For k = 
3, n + m - 1 with M2 = K - (M3 + ... + M„+m) we have 



(fc) 
ay = e 



r^n + m ^-*-fm + ---+*fn_|.m (m) 

(fc) , ^m=k + l ^i^l "i 

3 ^ A^fc + --- + *fn+m 



j = l,...,(Mfc + ... + Mn+J 



1 < < < ... < nZU,,,^M„^^ < M^-, + ... + M„+^. (4.44) 
Therefore, the function is the eigenvector of the initial permutation problem (14.271) 



M 

Qi^&MxeQPeeM i=i 
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where we denoted explicitely the 'empty sites' (with no excitation) as 0o(2 |). 

The coefficients ^^^\id,Q) gathered in ^^^\id) are connected with the next level coeffi- 
cients (id, Q). To see it, we fix the value of M3 and consider vectors B which characterizes 
excitations of the form 



M-Ms Ms 



5 = ((n + m) r,...,(n + m) T, 4 T, 4 T, 3 T, 3 T) • (4.46) 
Then, we consider the restriction of the relation (14.351) to 

^^'\zd) = V<l>(3)[^d,g,5]|i?,,,..,5,,, > . (4.47) 

restric. — _ 
Q,B 

One can recognize in this term the (B) coefficient (here we take only first M — M3 

values B): 

M-M'i 

5^<i.(=')Kg,5]efi,(,)...efi-;;i,,3)^ H <;0o(3t) = 0SU(5), 

(4.48) 

with Ui = q~^{i) for i = 1, M — M3. In the left hand side of the equation e^-i(j) are the 
operators which create the corresponding excitations Bi on the site q~^{i) of the chain of 
particles e^^. 

Therefore, using the same ansatz as in ( I4.30p for ^{y) in (j)^^_]^j,^{B) we can identify the 
coefficients ^^^\id,Q) as 

M-M3 

$(3)(zd,g,5)= '^'^'^P^Q'P^'^B) Hia^f, (4.49) 

with Q' G &M-M3 defined by q'{i) = q{i) for z = 1, M — M3 and 

5, = 4 t, (n + m) T, z = 1, M - M3 . (4.50) 

In the general case, the coefficients ^^^\P^Q) are defined by the same relations: for 
/c = 3, n + m — 2, we have 

Mfc+1 + ...+Mn+m 

$«(^rf,g,5)= <^^'^'\P,Q'P-\B) H (4.51) 

where Q' e 6^4+1+. ..+Mn+„ is defined by g'(i) = q{i) for i = 1, Mk+i + ... + M^+m, 

B, = {k + 1) t, (n + m) T, i = 1, Mk+i + ... + M^+m 

and there are relations similar to fl4.33p and (14.341) . 

At last, when A; = n + m — 1, using the results of gl{2\l) © gl{2) model, we get 

Mn+m-l+Mn + m 

<|.("+--^)(zci,g,5)= Y n [4"^]^- (4-52) 
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Equations of the type fl4.49p and fl4.5ip together with relations fl4.33p and (14.341) allow us to 
derive all the coefficients ^^^\P,Q). The eigenvalue reads 

i=l i=l 

The Bethe parameters n^'^), k > 3, ensure the correct multiplicity of eigenfunctions. Indeed, 
the total number of states 0m (^'•^\ n^^^ n^^'^"^^) is 

m+n— fc— 1 

V 1 = ; -, where Jk = Me, (4.54) 

, ^ (fe) (fc) ^ , 
— I 1+ 1 — ^ 

which shows that the ansatz is complete. 



5 Another gl{2\2) gl{2) model 
5.1 Comparison with AdS/CFT models 

It is known that the Hubbard model can be connected to the SU (2) subsector of the super- 
Yang-Mills (SYM) theory. The dilatation operator in this subsector can be identified with a 
Hamiltonian that is very close to the Hubbard one. Although not exactH, this correspondence 
has shed a new light on the integrability aspect of SYM models. In fact, one can introduce 
perturbatively a scattering matrix (obeying the Yang-Baxer equation) that differs from the 
Hubbard one by a phase (so-called 'wrapping problem'). This phase is also seen in the Bethe 
equation of the model. 

The models we have presented up to now possess the same property: they have Bethe 
equations that are very close to the Hubbard Bethe equations, but a phase. Unfortunately, 
this phase is built on 'hidden' Bethe parameters, but not on the impulsions of our 'particles' 
(which is the case of SYM Bethe equations). Hence, the correspondence is not immediate, 
but the present construction gives a way to introduce a phase in the equations. To strengthen 
the present approach, we show in this section an example of another 5f/(2|2) © gli^) model 
using different choices of the projectors tt and vf . It will lead to Bethe equation with a phase 
that (partially) depends on the impulsions of the particles. 

We recall that the Hamiltonian is given by 

L 

Hgl{2\2)®gl{2) = ^ ((SP)t x,x+l + (X'P)l x,x+l + U Ci^^ , (5.1) 

x=l 

where we choose the projectors Tr^ and 7r| such that A/| = {1, 2} and ATj = {3, 4}. Thus 

(SP)t .,.+1 = E E (^f^^tUi + ^f.^f.+i) ' (5-2) 

i=l j=3 

(SP)| ^^^+1 = Ef^Ef^^^ + Ef^Ef^^^ , (5.3) 

n — _L ??22 77133 77144 . ^ _ 7711 7722 (r A\ 

■^The correspondence breaks down at level 4 of perturbation theory [18] 
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We use the same approach, i.e. the coordinate Bethe ansatz, to solve this modeL Here 
we do not give all explicit details of calculation. The method is the same, however there are 
some modifications appearing when we pass from the initial problem to the first auxiliary 
problem and then to the second auxiliary problem. We briefiy give the most important 
statements as well as the Bethe equations. 



5.2 Bethe equations for gl{2\2) gl{2) 

As we shall see in the next subsection for construction of the Bethe ansatz, in this model 
there are four different kinds of "particles" above the vacuum state: e^^ which is defined by 
projectors above as vr-particle and e^-'^,e^^ and e^^^ which are the 7f-particles. To define the 
Bethe equations, we first introduce 

A = {ai, a2, . . . , a^i} for some integers such that 1 < ai < 02 < ... < a^-^ < N . (5.5) 

Then, the Bethe equations can be written as 

N2+N3 ■ . I, _|_ , u 

e"^-^ = (-1^^-^"^^"^^"^) n ]Zkl7 je[l,iV]\A (5.7) 

I bill rvj ~r f'^Ttl — 



1 rr ^sin/c, +z6^ + f i-r ^.ik, ""t^' - ih + % 



for m = 1, . . . , iV2 + 
k{n) = exp ( ' ^ -^1 < ••• <^7V3 <iV'2 + iV3 (5.9) 

where L is the number of sites considered in Hubbard model, N is total number of all 
e^^^e^^ and e^^^ "particles". Ni counts e^'^ excitations, N21N2, count respectively e^^ and 
e"^^ particles. Remark that with respect to the Bethe equations computed in the previous 
sections, the phase A(n) has been changed to 

A(n) ^ K{n)We-'^^ 

jeA 

showing a (partial) dependence on the momenta of the particles. This "dressing" of the 
phase is similar to the one suggested in [20]. 
The energy associated to the state is given by 

E = L-2{N - Ni) + 2 cos{kj) (5.10) 

je[i,N]\A 

and the momentum reads 

N 
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The set of integers A is related to the vr-particles in the first auxihary problem and the 
integers Uj correspond to the Bethe parameters of the last level, but their Bethe equations 
have already been solved: they just correspond to the quantization of these parameters. 
The parameters {&/}/=!,. ..,7V2+Ar3 = {0'Ni+i}i=i,...,N2+N3 correspond to 7f-particles in the first 
auxiliary problem, and do have Bethe equations, see eq. fl5.8p . 

5.3 Calculation description for gl{2\2) ® gl{2) 

In this paragraph we briefly describe some important points of the approach for this new 
model. At the first level of coordinate Bethe ansatz, together with the 7f-particles (the 
"physics" of which we studied above), we include some vr-particles. The states of N excita- 
tions can be written as 

m = Y.'^[^^'^Ki-4: (5.12) 

X 

withyl = (2,T); (3,T); (4,T); (2,i). 

This modifies the ansatz for the wave function as, for Xq^ < Xq^ < . . . < Xq^^ 

^J*(x)= HPQ,P~')e'^''"', P = P.^P' (5.13) 

with the energy: 

EP^. = 2 ^ cos(fcp^,(o) +L-2{N- iVi) (5.14) 

where Q G ©at. We have to consider the permutation of Bethe roots kj in some factorized 
form: P' = PttP^ where the terms permute only vr and vf particles separately. In addition, 
we could also vary the value of energy by mixing the impulsions of all particles, adding a 
permutation P^^ in the term e"^^^'^, but it does not produce new (independent) eigenvectors. 

Applying the Hamiltonian fIS.ip on the vector 05.121) we find the relations between the 
coefficients ^{PQ,P~^). Again we can gather all relations in a vector $(P) 

^P')^Y.^^{P',Q') |Ai,...,A^> 

Q',A 

where we have defined P' = PQ E &n, Q' = P~^ & &n and the sum is over all types of 
excitations and all corresponding sectors. The vector \Ai, ...,An- > belongs to Vi <^ ... ®Vn, 
where V = span{2 |,3 |,4 |,2 |} and represents one type of N excitations. 

Thus, the relations between ^{PQ, P~^) can be expressed using S'-matrix presentation 

^Ui2P) = s{'^{Xi-X2MP) 
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where >S'j2''(Ai — A2) is 



'S'i2'*(Ai — A2) - 
1 

ii2 



''12 



ri2 



ri2 



-iki 



ri2 



!'12 



-1 



(5.15) 



ri2 



'-12 



-1 / 



Expressions for ti2 and were given in previous sections, for example in fl3.26p . 

The periodic boundary conditions can be written using the S'-matrix and the vector $(P) 

^PC) = e^'^'^^^^P) (5.16) 

and if we choose P = C^~\ we arrive to the first auxihary problem. Thus, for j = 1, ...,N, 

Sf^ry..S^S'ff...S^,^^izd) = e^'^'^Hzd) (5.17) 



still with the convention S 



(1) 

ab 



The eigenvectors for this auxiliary problem are given by 



xe[l,Ar] 



(5.1^ 



with Ai = (2, 1); (3, t); (4, |) on the vacuum state filled by e"^^ particles. We recall that A''i 
counts e^^ particles and A'^2, ^3 correspondingly e^^ and e^^ particles. 

Comparing with the previous cases, the sector with two types of excitations e^^^ and e"^^ 
have been already treated but we have an additional vr-particle e^^. We write the eigenvector 
of this excitation similarly to (13.44^ in section 13. 4t 



N 



0[2 T] = E^= 



Aa e. 



2T 



3^^, with hx{a) = 5{x — a) 



(5.19) 



x=l 



This form of eigenfunction is supported by the fact that the eigenfunction in f l5.17p should be 
independent of index j. The ansatz for general case with for all types of excitations e^^, e^^ 
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and e^^ can be written as, for ^^(i) < Xq(2) < ■■■ < Xq(Ni+N2+N3) and Q G &N1+N2+N3 

Ni N2+N3 

*q(x,a)= 5^ $j«,n^-«(«") n (5-20) 

PeVf i=l n=l 

where Vf is the set of permutations acting on vf-particles only, and is defined in fl4.19p . 

The Bethe parameters ai, ...,aNi being the arguments of hx{a) are aheady quantized on 
the small chain [1, A^] and we choose them as: 

1 < ai < ... < oat-^ < iV. 

There are two different cases possible: 1) there exists a Bethe root Oq, = j for some 
a G [1, Ni], with j being the index in (15.1 7p and 2) there is no such Bethe root. In the first 
case, functions flS.lSp are eigenvectors of (15.171) with the eigenvalue 

sf:i,^...S^i]s[f...sf\^zd) = (-l)^-^i+ie*'=^(^^+^«)«l>(^ti) . (5.21) 

As we can see, the remaining parameters with i G [A''i + 1, A''i + + -^3] as well as 
the coefficients ^p'^i are not constrained in this first case: their Bethe equations is obtained 
from the second case. 

In the second case, when there is no Bethe root aa = j for any a G [1, A^i], (15.170 implies 
the conditions on the coefficients $p?i 

$p?l - ^«p(<?(i+l)) + f ' 

for all p(g(z)),p(g(z + 1)) in [A^i + 1, + + N3] and 



p-i —ik, 



e 



for all piq{t)) G [1, iVi], p{q{i + 1)) G [iVi + 1, iVi + iVa + iVs] and ai{a) is defined in (lOTD . 
The calculations for the eigenvalue of (15.181) give 



N2+N3 



S)%r--Si^]s\i^...S\^^^^^^ n ^A^^^N.)m). (5.24) 



1=1 



The periodic boundary condition on the coefficients is 

PQCJVJ+JV2+JV3 N 



'^)^'YlMap{q{Ni+N2+N3))) (5.25) 



1=1 

for Q such that q{Ni+N2 + N3) G [A''i + 1, Ni + N2 + N3]. CN1+N2+N3 is a cyclic permutation, 

Ck = nj,...nf-\ 

The difficulty in this model and of all models with vr-particles is that starting from the 
auxiliary problem we can not mix the Bethe roots of different types of excitations. Indeed, 
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in (15.201) the permutation acts only on vf-particles, while $(-P) mixes a priori any kind of 
particle. In this case we define again the vector $(-P), but only for vf-particles as actually it 
was defined in fl3.72p : 

<I>(P)= $J(A)|A,„...,A,(^,+;v3)> (5.26) 

where the sum is over all types of vf excitations and all corresponding sectors. The vector 
\Ag^, Agj^ > is in Vi ® . . . (g) V/^ with V = span{e^^ e'^^}. 

Working with the periodicity condition fl5.25p . we take for instance Q = id, which ensures 
that the constraint q{Ni + N2 + N3) G [iVi + l,Ni + N2 + N3] is satisfied. Thus, in the left 

hand side, in the coefficient ^^^l^'^'^i+^a+'^s) p(n^+n2+n^) identify which particles are 

permuted by n^|^^_|_^^_|_^^^ and we can apply corresponding relation (15.221) or (15.231) . Due to 

the fact that P G "P/ does not act on vr-particles, we can use relation (I5.23p to extract from 

■p{j) 



the coefficient the action of the last Ni permutations n!'|"'2 :M_i_M\,j = ^,---,Ni, and get: 



iVl T-,p(Ni + l) Y,P(Ni+N2 + N3-l) N 

i=l 1=1 

(5.27) 

We can rewrite this condition using $(P) notation. Choosing P as a power of the cyclic 
permutation acting only on 7f-particles, 

p _ /-iN2+N3-m _ YtNi+N2+N3-l\N2+Ns,-m 
^ — ^ — \^^Ni+N2+N-i---^^Ni+N2+N-i ) ' 



the calculation becomes equivalent to the one of section 13.41 and we obtain: 

Ni N 

JIe-^^-"«a„,(a„)<l>(C^^+^^-™+i) = (-1)^^ J] aKa^)$(C^^^-^^-'") • (5.28 

i=l 1=1 

Thus, for m = A^i + 1, iVi + iV2 + A^3 

Me* °'0"ai(am) I I I -. : — I {Pmm.+l---Prn Nx+N2+N3Pm.l---Prnrn-l)^{'id) 



N 



1=1 

and finally we get 



^{id) = 0, 



e^-i^n--"- n ( ' '"^! )-(-i)"-n-^'(^- 

i=l i=JVi+l \ ' ™ ' 2 / i=l 

m = A^i + l,...,A^i + iV2 + iV3 (5.30) 
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where we have introduced the set A = {ai, 02, ... , aN^} C [1, A^] and the Bethe parameters 
rij that label the eigenf unctions of the cyclic permutation (as in section 

To stress the difference between the quantized parameters aj, j < Ni, (that are similar 
to the parameters Uj), and the parameters aj, j > Ni, we denoted the latter bj = clj+Ni in 
the Bethe equations written in section 15.21 



6 gl{2\l) e gl(2) and gl{2\2) gl(2) Hamiltonians 

In previous sections we considered examples of Hubbard model with algebra gl{2\l)^ ® gli'i) i 
and its generalization to gl{x\\m)^ (B gl{2)i model. However only examples with "small" alge- 
bras like gl{2\l)-^(Bgl{2) I or 5f/(2|2)^ ©(yf/(2)j^ seem to find applications in physics. Performing 
different Jordan- Wigner transformations one can write the corresponding Hamiltonians in 
terms of creation and annihilation operators. 



6.1 Jordan— Wigner transformation 

We briefly recall some relations of Jordan- Wigner transformation [21] (for more detailed 
explanations, see e.g. [16]). The Jordan- Wigner transformation essentially consists in the 
construction of a mapping 

E'^ e gl{2P-'\2P^') ^ {c\ c ; d\d; e^ e ; ...}, 

where c,d,e,... are fermionic operators. To present this construction, it is convenient to 
introduce a matrix X 

^=(1/ n= = ctc. (6.1) 

Its entries X^^^, (/i, z/ = 1, 2) have a natural grading given by [fi] + [u] where [1] = 1 and 
[2] = 0. The mapping is given by the relation 

E'^ - ® 41 ® - ® 41 (6-2) 

where to every value i and j G [1,2^] is associated with an element {fii, /i2, ... ,/ip} 
and {z/i, z/2, ... , t'p} respectively with /Xj, i^j = 1 or 2. Total grading is given by s = 

Transformation for gl{2)i algebra is simply given by identification of matrices E'''^ and 
Xfj^i^. But since gl{2) algebra contain only bosonic operators, in order to satisfy the an- 
ticommutation relation between fermionic operators on different sites {{cl,Cb} = 0), one 
should introduce some factor to E'^J (elementary matrix i?*-' on the site a). To present the 
result, we gather the different matrices i,j = 1,2 into a formal matrix E. Then, the 
Jordan-Wigner transformation for gl{2) algebra reads 

E„=(|; S)^C7'' n°') n (6.3) 

Now consider an example for the cases of gl{2\2\ and gl{2\V)^ algebras. For both algebra 
we perform one mapping and in the case of smaller algebra {^gl{2\V)^ we remove a subspace. 
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We take the mapping {1 ^11, 2 ^ 22, 3 ^ 12, 4 ^ 21} as an example. In the initial 
problem, elements oi gl{2\2)-^ (B gl{2) i and gl{2\l)^ (B gl{2)i algebras transform into fermionic 
operators which are considered to anticommute even for different spins, e.g. {c|,C|} = 0. 
This implies that the total Jordan-Wigner transformation is written as 



(1 - KMit nt cn 



c d 



\ 



n 



;i - <^)4t 



n(l-2n 



(6.4) 



a;=l 



(6.5) 



forE=(E*^),,,=i,„„4e^/(2|2). 

For E = (-£'*■' )jj=i, 2,3 G fi'^(2|l), it is given by the same matrix without line and row 3 

/(I -<t)(1 -^at) -(^a^da^ -Ca^{l - n'i^)\ L 

EaT = cl^dl^ <T^aT '^aT^T 11(1 - 2<J 

where = 6j, „ fe^^ ^ is the particle number operator for b = c,d, and a =T, i- We have also 
standard relations between the operators: 

X ^ C/s y} = S^ry ] {dl ^ , y} = S:,y 6^13 ] {Cax,dpy} = 0, J.. (6.6) 

Note that the operators d^ are not present in the construction, so that one can drop the 

arrow on the operators dp d^. 

One can remark that the different choices of mappings on {/i, z/} is equivalent to 
some transformations on fermionic operators' level (e.g. — > c, etc...), therefore all the 
Hamiltonians are equivalent in this sense and differs one from another by changing the 
representation (e.g. from electrons to holes). 



6.2 gl{2\2) e gl{2) model. 

The Hamiltonian of the model is given by 

L 

Hgl{2\2)®gl{2) ■ 



with 



Huuh + ^{c\ ^+iCt ^ + c{ ^ct ^+i) (nf-^nf +^ - ^ - +^) 

x=l 

L 
x=l 

L 

-2n^(l-2<^)(l-<^)nf^ 

L 

+ x+i'^T x + d\ x+i) ('^xt''^x+it ~ ""-xT ~ '"'x+it) 



a;=l 



X=l 



(6.7) 



(6.^ 
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The eigenfunctions for this Hamiltonian are made of creator operators c|, and c|. They 
can be written in the following form and correspond to the solutions found in the previous 
sections: 

Ni-n Ni+N2-n n N3 

<u,^3 = EE E n-^^^niit n <Tn4T<Tniiio> 

? M 2 j = l j=Ni + l k=l j=l 

-k^H Vk^vi ^i^^j'Vk 

= , with n = 0,...,min(iVi,iV2) (6.9) 

and A = (4Tr^,3T^, 2X^^,277^?^ (6.10) 

<f)[A] is the eigenfunction given in fl4.10p . Remark that, in addition to the particles c|, (i| 

and c| (corresponding to 4 t, 3 t and 2 | resp.), we have a doublet c|(i| corresponding to 

2 f. The particles cj and c| can be identified with a (spin up and down) electron, while (i| 
can be viewed as a spin fermion that can form bound state with the spin up electron. 
The energy of the excited state ^j^^-* reads 

Ni+N2+N:i-n 

= L- 2{Ni + N2 + N3-n) + 2 ^ cos fc, , 

1=1 

where the parameters ki are Bethe roots defined by equations given in section 14.11 (with 
n = m = 2). 

6.3 gl{2\l) e gl{2) model. 

The Hamiltonian of the model is given by 

L 



HgH2\i)®gi(2) = + E (4^+1 ^ + 4^^T ) (^T x^x+iT -'^Tx -^T^+i) 

x=l 

L 

+ E (4 x+l^l xd\ 3;+iC?T a; + 4 x^T x+ld\ x+1^ 
x=l 
L 

-«J](l-2<.^)(l-0?, (6.11) 

x=l 

where Huub has been given in (16.81) . Again, the eigenfunctions for this Hamiltonian corre- 
spond to the solutions found in the previous sections. They have the form (for N2 < Ni) 

Ni N2 N3 

^N.,N2,N3 = E E E ^(^'^"''^"')niTniT<Tn4,iio> 

S y X 1 = 1 1 = 1 1 = 1 

^k^H Vk^yi ^iT^^j'Vk 

Ni N2 N3 

= (-l)^^+^^cp[A], with A = (3Tr^,2Tr^,2Tr?^) (6.12) 
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(f)[A] has been defined in f l3.14p and the corresponding eigenvalue reads 

N1+N2+N3 

E = L- 2{Ni + iVs + iVs) + 2 ^ cos k 

1=1 

with Bethe roots ki obeying the equations given in section I3.2[ 

7 Conclusion 

In this paper we presented the Bethe equations using the coordinate Bethe ansatz for gl{2\l)(B 
gl{2) and gl{y\\xa) ®gl{2) generalized Hubbard model. We wrote explicitly the Hamiltonians 
for several cases in terms of fermionic creation and annhilation operators. Clearly, the full 
derivation of the Bethe Ansatz Equations for generalized Hubbard models [16, 17] has to 
be accomplished: the case of 5f/(n|m) © (7/(3) and its generalizations to gl{x\\m) © gl{x\!\x(\.') 
algebras are presently under investigation. 

Applications to condensed matter physics deserve also some attention. The models pre- 
sented in section [6] are generalization of Hubbard models to several types of fermions. They 
could be of some relevance to systems where electrons of different 'types' or 'colors' occur 
(for instance on some ladder spin chain). 

Although the link with AdS/CFT correspondence is not direct, the present construction 
gives a way to introduce a phase in the Bethe equations of Hubbard type models (see 
discussion in section ISTTI) . It is thus worthwhile to look deeper at these models. They could 
be good candidates for an integrable model close to the one underlying the SYM theory, or 
give a new point of view for the wrapping problem. 

Obviously, in all cases, the thermodynamical limit of these models needs also to be 
investigated. 
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